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Abstract 

We introduce a large class of scalar-tensor theories where gravity be- 
comes stronger at large distances via the exchange of a scalar that mixes 
with the graviton. At small distances, i.e. large curvature, the scalar is 
screened via an analog of the Vainshtein mechanism of massive gravity. 
The crossover distance between the two regimes can be made cosmolog- 
ical by an appropriate choice of the parameters. 
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There has been recently a renewal of interest for theories where gravity is modi- 
fied at large, cosmological, distances. This is due in part to the wish to find expla- 
nations to the observed acceleration of the expansion of the Universe not relying on 
a new dark component. One prototypical example of such theories is DGP gravity 
[Tj with its interesting cosmological consequences |2| 13|. Other motivations exist, 
like the dark matter puzzle |1] or recent observations made in Ref. [5]. It is how- 
ever notoriously difficult to modify General Relativity (GR in the following) at large 
distance in a consistent way without spoiling the very good agreement between GR 
and various experiments and observations, in particular those of the motion of solar 
system bodies as well as strong field tests relying on binary pulsars (see e.g. |6j). 
Besides DGP gravity (which suffers in particular from unsolved questions related 
to its UV completion), degravitation models and some related proposals [7], some 
other classes of models have been considered with some details: e.g. the Chameleon 
models [8], or the - equivalent to scalar-tensor - f{R) [9]. What we really need is 
some theory that behaves as GR in the strong field limit as well as for weak fields 
appropriate to describe gravity in the solar system, while it deviates from GR in 
some ultra- weak field limit. 

A theory which was thought to have similar properties is nonlinear massive 
gravity, first considered in the context of strong interactions ^0\. It can be defined 
as a bimetric theory, where one of the metrics, say g^u, is dynamical (with an 
Einstein-Hilbert action) and minimally coupled to matter, while the other, say /^jy, 
is non dynamical and couples to g^^ (see e.g. [H]). The couphng between the 
two metrics is chosen such that, expanding around fiat space-time, one recovers at 
quadratic order the Pauli-Fierz action |T2], the only consistent action for a massive 
spin two. A massive graviton has 5 propagating polarizations, among which a scalar 
mode responsible for an extra attraction felt by non relativistic sources. This scalar 
mode, if one only uses the quadratic Pauli-Fierz action, leads to order one difference 
in PPN parameters from those of General Relativity, irrespectively of the smallness 
of the graviton mass (this is the famous vDVZ discontinuity p3]). However it was 
argued some time ago that the scalar polarization could hide itself in strong enough 
field configurations via the so-called Vainshtein mechanism [H]. This relies on the 
observation that the theory has hidden derivative self interactions in the scalar sector 
that shuts off the effect of the scalar attraction at distances smaller than the so-called 
Vainshtein radius Ry [HI [IE]- This radius, defined by Ry = m~'^/^R^J^ in terms of 
the graviton mass m and the standard Schwarzschild radius of the source Rs, can 
be very large; e.g. it is much larger than the solar system for a graviton of Hubble 
radius Compton length. The Vainshtein mechanism was argued to fail in nonlinear 
massive gravity |I7l[IH], a theory which anyway suffers from ghost-like instabilities 
[H]. However, it was recently observed exphcitly [|20], building on some previous 
qualitative arguments [H], ED, [22], that it does work in some limit of this theory 
where one only keeps the dominant derivative self interactions of the scalar sector, 
the so called "Decoupling Limit" (DL in the following). 

In this essay, we propose a general structure for a large class of scalar-tensor the- 
ories, in which the scalar field "camoufiages" in strong enough gravitational fields. 
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via a derivative self-interactioiij^ Note however that to obtain this effect, the grav- 
itational field does not need to be very strong and, similarly to massive gravity, it 
can already happen in the solar system. The structure of our scalar tensor theories 
follows closely the one of the DL of nonlinear massive gravity and, before introducing 
the former, we remind some crucial properties of the latter. 
In nonlinear massive gravity, the metric field equations read 

G^, = (T^, + T^g , (1) 

where G^,^ denotes the Einstein tensor computed with the metric g^^^, T^^, is the 
matter energy momentum tensor, and T^^, is an effective energy momentum tensor 
coming from the interaction between metrics. It depends non derivatively on /^,^ 
and g^y. Taking a (7— covariant derivative of Eq.Q we obtain the constraint 

V^Tf^y = (2) 

which T^^ should obey. Following Ref. [HI [TT] one can then look at spherically 
symmetric solutions, aimed for example at describing the metric around stars, using 
the ansatz 

g^^^dx^'dx" = -e^^^^dt^ + e^^^^dR^ + R^d^^ , (3) 
f^ydx'^dx'' = |^i_M^ye-^(^)cii?2^e-^(^)i?W, (4) 

where a prime denotes a derivation w.r.t. the radial coordinate R. With such 
an ansatz, g^^, is easy to compare with the standard form of Schwarzschild solu- 
tion, while ffj,i, describes a Minkowski space-time in some unsual coordinate system, 
parametrized by the function fi to be determined. In the DL, the system of equation 
to be solved, Eqs. ([Tp| collapses to [20] 



^ + ^ = -^^'(3/i + V) + ^, (5) 



R~R^ 
R^~ 2R 



m^H, (6) 

Qip), (7) 

+p + 5Rfi'^ + ^Rfifi" + ^R^fi'fi"^ I , (8) 

where p represents the source energy density, and Q{fi) contains the only left over 
nonlinearities parametrized by a and p. This parametrization covers all the pos- 
sible interaction terms between the metrics [20j. In the DL, those nonlinearities 



^Hence the title of this essay, where "k" refers to "kinetic". 
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only appear in equation ([2]), which is there equivalent to Eq. ([7|. The first two 
equations ([5]|6]) correspond to taking the DL in equation ([T]). corresponds to 

the strongest scalar derivative self interactions of the model |20]. These interac- 
tions, when expressed in term of canonically normalized fields, are suppressed by 
a scale A = m 

4/5^1/5 

|16| that can be made here explicit by rescaling the fields 
as z/ — > Mp^u, X Mp^X, fx — > Mp^m~'^^, in which case the right hand side of 
equation ([T]) simply appears as divided by A. In fact the DL is simply obtained 
by using those rescaled fields in the original equations ([l]j2| and letting Mp oo, 
m — > 0, while keeping A (as well as T^^/Mp) fixed. It is expected to give a good 
description of the solution in the range Rq <^ R <^ where Rq is a distance 

scale that can be made parametrically lower than Ry, and even extend down to the 
Schwarzschild radius [221 [20] • The structure of equations ([5]j7]) can be understood 
from the fiat space-time action 

f r 

S=^ d^x^2h'"'d^d^h - 2W''dyd,h1 + W^Uh^^ - hUh 

+ [AQi^.d^'d'ct) - hUcj)) + 4a (□0)^ + 4/3 (00 0,^, (j)^^')] } (9) 

for the dynamical field h^y and 0. Indeed, using the identifications [20] 

V = {A, 4, (10) 
/i = -207i?, (11) 



h^u field equations correspond to Eqs. while the one for correspond to Eq. 

([7|. The peculiarity of action ^ is that does only get a kinetic term via a mixing 
with h^y [H], this being entirely due to the structure of the Pauli-Fierz mass term. 
Outside of a source and at large distances, the nonlinearities Q(/i) can be neglected, 
and one finds that the system ((5j(7| is solved by 



C C I C 

where C is a constant of integration expected to be proportional to Gtv and to be 
fixed by matching to the source. Notice that one has v ~ — 2A, in contrast with the 



GR result, this being due to the scalar exchange. The behaviour (12) is vahd for 
R ^ Ry- Below Rv, inspection of equations ([5]j7| shows that the nonlinear term 
Q{fi) dominates over the terms linear in /i and at the same time is of the order of 
X'/R ~ X/R^ ~ v'/R. Hence, for R <C Ry, one can neglect the terms hnear in 
in front of those in the left hand side of equations ([5| and ([6j, and as a result, one 
recovers linearized General Relativity (obtained for = 0). This is the essence of 
the Vainshtein mechanism, and it was explictly shown to work this way by solving 
numerically the system of equations ([5]j7| [20]- Having this in mind, one can design 
a large class of scalar tensor theories in which one expects to have large distance 
modification of GR, but a small distance recovery a la Vainshtein. 
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Indeed, a covariantization of action ^ is given by 

S = Mj. [ d^xy/^ ( ^ + Jm^^i? + m'^H{(P) ) + S^, 



fl3l 



where H^cp) is some derivative (covariant and higher than quadratic) self interactions 
for that will be discussed below, and Sm represents the action for matter assumed 
to be minimally coupled to g^j,^ (7 is some order one parameter). These theories 
can of course always be written in an Einstein frame. When expanded around flat 
space time and = 0, action (13) reproduces ^ at quadratic order. It also contains 
derivative self interactions of similar to those left over in the DL of massive gravity 
and encoded into the function H. Some possible choices for H are following: 



H{(f))MG 

H{(j))DGP 
H{(t>)K 



|(n03) + ^(nv.v^;^., 

m2n00.^0''^, 

K{X), with X = m^(j).^(j)'>' 



rr? (0;A 



) [2 (D^ 



CovGal 



m [(p,x> 



2(00) -2 



R 



(14) 

(15) 
(16) 
(17) 

(18) 



The first choice is reproducing the structure of the self derivative of arising in the 
DL of massive gravity; the second choice does the same for DGP gravity [231 [2lj5 
the third possibility, Hk (chosen such that dxK{X)\x=o = 0), corresponds to K- 
essence-like model^ [25]; the fourth case corresponds to the Galfleon model [26] , 
while the last case is the same model with an appropriate non minimal coupling to 
curvature removing higher derivatives in the e.o.m. [27] (see also [28]). The metric 



field equations obtained from (13) read 



G^^ (1 + 7m V) + im^ {g, 



(19) 



where Tj^^ comes from varying H{(f)). The scalar field e.o.m. is 



;7R + £4 



(20) 



where £^<^ come from varying H{(j)) w.r.t. 0. What is first interesting to notice is 



that aU the models defined by (13) and their corresponding H have a "Vainshtein 



radius" and a "decouphng limit" which isolates the strongest scalar self interaction. 
This interaction reads schematically 



N 



f 



^Note added after this work was submitted to the 2009 Gravity Research Foundation Essay 
Competition (March 31st 2009): the recent arXiv submission |32] considers aspects of the cosmol- 



ogy associated to a model close to (15 I (see also [33j). 

^With however the important difference that in the Einstein frame, matter couples to a metric 
which depends on (j>. 
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with A'^^ and Nf integers. Using then the scahng (12) which also apphes here, as 
well as (10 11), one has 



H{^) fRv,H 



Nf+Na-i 

hD^ \ R~) 

where we obtained the Vainshtein radius corresponding to the model H as 



In the DL, and for spherical symmetry, equations (19) reduce to ipmn, while Eq. 



(20), once integrated once, leads to an equation analogous to lujr\ The DL can 
be used as a starting point for a numerical integration of the full nonlinear e.o.m. 
However, it is highly non trivial that non singular solutions found in the DL will cary 
into non singular solutions of the full theory, as we learned from nonlinear massive 
gravity |T7[I^. 

One of the main result reported here is that we were able to find non singular 



numerical solutions of the full system (19 20) with sources, for some particular 



choices of H. Those solutions show a recovery d la Vainshtein of GR below the 
respective Vainshtein radn. Our expectation is that such solutions will also exist for 
a large subclass of "k-Moufiage" type of models to be determinecQ Figure [T] shows 
for example some of the results of the integration of the full nonlinear system for a 
simple Hk given by Hk{X) oc as well as for Hdgp- Notice in particular that 
the HoGP model has the interesting property that it has the same Vainshtein radius 
as the original DGP model as well as PPN corrections inside the solar system. 

Various issues are left for future works. One obvious question, not to mention 
possible troubles with superluminal propagations (see e.g. the different points of 
Refs. |30l[3l]), has to do with the stability of the different model considered. E.g. 



it is clear that the choices (14) and (17) leads to higher derivative field equations 



(see [27] for what concerns (17)), and hence unstable modes. Higher derivative, 
however, would not appear for the other choices, but this by itself does not ensure 
stability of the model. It is also clear that the model considered here should be 
UV completed in some appropriate way. The situation is the same as for DGP 
gravity, Galileons, or generic models of K-essence and one should worry about the 
low scale associated with the dominant scalar interaction. We believe however that 
the framework presented here, if only phenomenological, should allow various novel 
effects of large distance modification of gravity to be investigated in a cosmological 
context, a task that is certainly worth pursuing. 



^The required integration comes from the relation (111. 

^Note that such non singular (and in this case higEer dimensional) solutions are not known 
explicitly in the original DGP model, where only approximate or DL solutions have been obtained 

CSIEIIMI- 
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Figure 1: Left figure shows the derivative of the scalar field, 0', in the presence of a 
source of radius R = 0.01 (in the unit of the Vainshtein radii of the models) for H = 
Hk{X) oc (dashed blue curve), and H = Hdgp (red solid line). The different 
asymptotic regimes are shown by thin fines. One sees in particular the transition 
happening at the Vainshtein radii Rv,h- The right plot shows the derivative of 
the functions A and z/. A' (dashed blue curve) and z/' (red solid curve), appearing 
in the spherically symmetric ansatz ([3j, along with 0' (thin black curve), for the 
H = Hk{X) oc model and outside the source. One sees the transition between 
GR regime {R -C Rv,h) and the scalar-tensor regime {R 3> Rv,h)- Below Rv,h the 
scalar field "camoufiages", i.e. its contribution becomes subdominant. 
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